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Model analysis on thermal UV-cutoff effects on the critical boundary in hot QCD 
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We study the critical boundary on the quark-mass plane associated with the chiral phase transition 
in QCD at finite temperature. We point out that the critical boundaries obtained from the Lattice 
QCD simulation and the chiral effective model are significantly different; in the (Polyakov-loop 
coupled) Nambu-Jona-Lasinio (NJL) model we find that the critical mass is about one order of 
magnitude smaller than the value reported in the Lattice QCD case. It is known in the Lattice 
QCD study that the critical mass goes smaller in the continuum limit along the temporal direction. 
To investigate the temporal UV-cutoff effects quantitatively we consider the (P)NJL model with 
only finite A'^,- Matsubara frequencies taken in the summation. We confirm that the critical mass 
in such a UV-tamed NJL model becomes larger with decreasing Nt, which demonstrates a correct 
tendency to explain the difference between the Lattice QCD and chiral model results. 
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Introduction The phase diagram of quark matter as 
described by Quantum Chromodynamics (QCD) is of 
great interest in theoretical and experimental physics. 
Since QCD has been established as the fundamental the- 
ory of quarks and gluons, it is our ultimate goal to un- 
derstand all the phenomena related to the strong interac- 
tions from the first principle theory. It is difficult, how- 
ever, to surpass perturbative calculations due to largeness 
of the QCD coupling constant in the low-energy regime. 
One then has to make a choice among non-perturbative 
techniques such as the Lattice QCD (LQCD) simulation, 
some effective model descriptions sharing the same global 
symmetry as QCD, etc. The Lattice QCD (LQCD) is a 
theoretical framework of quarks and gluons in discretized 
space-time JL] which is suitable for non-perturbative in- 
vestigations. The Monte-Carlo simulations in the LQCD 
have been recently developing to a reliable level not only 
for hadron spectrums at zero temperature but also for 
thermal properties in hot QCD matter 0. 

Critical phenomena in hot QCD can be studied in 
the LQCD simulation [3|] and it is an important issue 
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to determine the order of the chiral phase transition as 
a function of the temperature T and the quark chem- 
ical potential fi. In this paper, we study the critical 
boundary in the chiral phase transition on the quark- 
mass plane, as represented by the "Columbia plot" [1]. 
Figure [T] shows a schematic picture of the Columbia plot 
at fi — where each point in the parameter space of cur- 
rent quark masses {rriu, rud, rris} is marked by its order 
of the phase transition when T is raised. Usually one 
sets niu = md{= niud) and studies the order of the phase 
transition on the mud-nig plane. Symmetry arguments 
show that it should be a first order phase transition in 
the chiral limit {iriud = m-s = 0) [5|, while for inter- 
mediate quark masses at which both chiral and center 
symmetries are badly broken, only a smooth crossover 
is possible. This defines the second-order critical phase 
boundary that separates the first-order and crossover re- 
gions. 

The "phys. point" in Fig. [1] designates the "real world" 
with physical current quark masses {m-ud = 5-5 MeV 
and TOs = 135.7 MeV in a model study), whose location 
should crucially decide the order of the phase transition 
at /i = 0. More importantly, the distance between the 
second-order critical boundary and the physical point on 
the Columbia plot gives a rough measure of how far the 
QCD critical point is separated from the fj, — situa- 
tion. Roughly speaking, the farther the physical point is 
located away from the critical boundary, the larger fi is 
necessary to make the critical surface hit on the physical 
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FIG. 1: Columbia plot at finite T and jj. — 0. 



point. 

The critical boundary can also be analyzed by means 
of the LQCD simulation since the Columbia plot is 
drawn at zero baryon density, where the Binder cumu- 
lant is a useful quantity to identify the order of the 
phase transition [a, 0, HI- To draw the Columbia plot 
we can use the three-flavor Nambu-Jona-Lasinio (NJL) 
model which is a model of quarks with the effective 
Lagrangian so constructed to have the same global sym- 
metry as QCD [l^jliii- (For other chiral model studies 
see Gauge symmetry is, however, absent due to 

the absence of gluons in the model. It is a crucial aspect 
in the NJL model to incorporate chiral symmetry as an 
important ingredient in understanding the properties of 
light hadrons. Further, the model can successfully de- 
scribe the spontaneous breaking of chiral symmetry at in 
the QCD vacuum and its restoration at high tempera- 
ture and/or density. The three-flavor NJL model gives a 
first order phase transition in the chiral limit by virtue 
of the Ua(1) symmetry breaking term while, at in- 
termediate quark mass, the chiral restoration results in 
a crossover due to explicit chiral symmetry breaking by 
finite quark mass as long as /i ~ 0. Thus, the NJL model 
should be able to reproduce the expected structure of the 
Columbia plot (i.e. critical boundary in the left-bottom 
corner in Fig.[T]), which can be compared with the LQCD 
results. 

Although the NJL model successfully describes the 
qualitative feature of the chiral phase transition, there re- 
mains a significant difference between the critical bound- 
aries from the NJL model and the LQCD simulations. 
In fact, the value of the critical quark mass is about one 
order of magnitude different, which is crucial in the QCD 
critical point search. One may have seen a scatter plot 
of the critical point locations from the LQCD atternpts 
as well as the chiral model approaches as shown in [l4| . 
However, there is a huge discrepancy in the critical mass 
boundary already at zero baryon density depending on 
the adopted methods, and then, how can one trust any of 
them at finite density? Let us pick some concrete num- 



bers up. We define the critical quark mass rric as a point 
where the niud = f^s symmetric line intersects with the 
critical boundary. Then, rric is around 14 MeV in the 
LQCD simulation [6*1 but it is only 1.1 MeV from the 
NJL model. 

In the LQCD studies 0, IBi '^^^ argued that the 
first-order region drastically shrinks by removal of the 
discretization errors. In [Sj, for example, the authors 
used the Symanzik improved gauge and stout improved 
fermionic action and concluded that the critical mass 
goes below 7% (and 12%) of the physical quark mass on 
the Nt = 4 (and Nt = 6 respectively) lattice. The con- 
clusion in is consistent with this qualitatively though 
the quantitative change is not such drastic. Because the 
number of spatial points is large enough (16, 24, etc) 
in the LQCD simulations, smallness of Nr in the tem- 
poral direction is crucial. In other words, the temporal 
UV cutoff ttNt-T has drastic effects on the critical mass 
boundary even at zero density (and hence on the QCD 
critical point at finite density as well, which is beyond 
our current scope in this paper). 

Here, we shall examine the chiral critical boundary by 
imposing a finite UV cutoff in the temporal (thermal) di- 
rection in the NJL model to mimic the LQCD situation. 
In this way, we can give a quantitative prediction about 
how much the critical boundary moves by the finite A^,- 
effects. In terms of Fourier transformed variables the 
summation over finite-A^T- temporal sites amounts to 
the Matsubara summation over finite-A"r frequencies. 
Hereafter, throughout this paper, we will drop the sub- 
script and simply write as A^ for notational simplicity. 
In what follows, we present a brief description of this 
modified NJL model and discuss the machinery of how 
to deal with finite A^ in this model treatment. Later, we 
shall also present some results using the Polyakov-loop 
coupled NJL (PNJL) model i, 0, [M Ezl, 13, for 
the PNJL model is a much better description of QCD 
thermodynamics near the critical temperature. 



Modified NJL model 
grangian is. 



The three-fiavor NJL model La- 



>Cnjl 9 (i^ - m) q + C4 + jC-e, 
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-gu [detqi{l - 75)9^ + h.c. ] . 



(1) 
(2) 
(3) 



The current quark mass matrix to in the kinetic term 
takes the diagonal form to = {mu-,rnd,ms), and we set 
TO„ = TTLd (= rriud)- In the above C4 expresses the 
four-fermion contact interaction with the coupling con- 
stant gs, where Aa is the Gell-Mann matrix in flavor 
space, and is the Kobayashi-Maskawa-'t Hooft inter- 
action [1^ [23, [2l| with the coupling strength gu ■ The de- 
terminant runs over the flavor indices, which leads to six- 
fermion interaction terms in the three-flavor case. This 
interaction explicitly breaks the Ua(1) symmetry, whose 
microscopic origin comes from instantons [2l| . 
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We should solve the gap equations which we can get 
by differentiating the thermodynamic potential with 
respect to the constituent quark masses; 



dm* 



= 0, 



dm*, 



= 0. 



(4) 



Here m* and m* are the constituent masses for u and s 
quarks. The thermodynamic potential is defined hy n — 
— lnZ/(pV), with the partition function Z, the inverse 
temperature /? = 1/T, and the volume of the system V. 

In the mean-field approximation for the NJL La- 
grangian ([T]) we obtain the following gap equations, 



rn* = ms + 2i gs NctrS 



2i5s N.trS" ~ 25DiV2(tr5")(tr5^), 



Here A^c = 3 is the number of colors and tr5* is the chiral 
condensate given explicitly by 



i trS' = 4m' 



(2^ 



(iT) 



n— — o 



E: 



■2 ' 



(6) 



where the Matsubara frequency is a;„ 



(2n 



1)ttT and 
A detailed 
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the quasi-particle energy is -Ei = -^/p-^ + 
derivation of the gap equations is clearly given in [1 

To study the effect of the finite cutoff in the temporal 
direction on the critical boundary, we then perform the 
following modification in the Matsubara frequency sum; 



E 



N 

E ■ 

n=-JV+l 



(7) 



We shall reveal the critical boundary in the mud-'m-s plane 
for various values of N. To this end, we adjust the model 
parameters to fit the physical quantities for a given N 
(each parameter set representing a modified NJL model 
with a different temporal cutoff). This is reminiscent of 
the familiar Wilsonian renormalization group (RG) flow. 
After the high frequency modes are integrated out, in 
principle, the coefficients of the operators run on the RG 
flow. Here we remark that we neglect interaction oper- 
ators in other channels than those in Eqs. ((T|)-(I3])) since 
their effects are small in the mean-field calculations per- 
formed below (see e.g., [10!] and refs. therein). 

The following are the model parameters in the three- 
flavor NJL model; 

current quark masses m^d, m^ , 

three-momentum cutoff A , 
four-point coupling constant i 

six-point coupling constant go ■ 

As for the light current quark masses, we fix m„d = 5.5 
MeV following [ii| , which gives isospin symmetric results 
and thus m* = m^. The remaining parameters are de- 
termined by a fit to the four physical quantities: the pion 
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rus (MeV) 


A (MeV) 




5dA^ 


15 


134.7 


631.4 


4.16 


12.51 


20 


135.0 


631.4 


4.02 


11.56 


50 


135.3 


631.4 


3.82 


10.14 


100 


135.4 


631.4 


3.75 


9.69 


CXD 


135.7 


631.4 


3.67 


9.29 



TABLE I: Fitted parameters for various values of A''. 



mass m^r, the pion decay constant the kaon mass mx, 
and the 77' mass m^/ , whose empirical values are 



m^ = 138 MeV , 
mK = 495.7 MeV, 



= 93 MeV , 
TO^' = 958 MeV . (8) 



Numerical results The limit of A'^ = 00 invariably 
corresponds to the standard framework of the NJL model 
and we then use the same parameter set fixed in [ii| . For 
finite N, we have picked up a temperature T = 50 MeV 
at which we have carried out our fitting procedure. This 
choice is just because of technical convenience. Besides, 
calculations at T = 50 MeV should be sufficiently close 
to those at zero temperature since the constituent quark 
mass is then several times larger than the temperature. 
The fitted parameters in this way for N = 15, 20, 50, 100 
at T = 50 MeV are displayed in Tab. HI 
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FIG. 2: Constituent quark masses mj^ (lower) and m* (upper) 
as a function of T for various values of N . 

As seen from the table, the couplings gs and gj^, be- 
come larger with decreasing N . If we fixed the coupling 
constants, the effect of smaller N would be to reduce 
the constituent quark masses. Therefore, to keep m* un- 
changed, we need larger 173 and to compensate for 
this reduction (see Eq. Note that, as expected, the 
whole behavior of m* , m* as a function of T monotoni- 
cally approach the standard (i.e. N = 00) case with in- 
creasing N . In Fig. [21 we only show the numerical results 
above T = 40 MeV because the small-T region is badly 
affected by finite N and the calculations are no longer 
reliable there. We can understand this intuitively; 1/T is 
the extent along the imaginary-time (thermal) direction 
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and thus, the extent becomes longer as T goes smaller. 
Hence, at lower T, reliable estimates need more N . This 
means, fixing the value of N , the choice of the sniall-T 
regions would automatically lead to unphysical results. 
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FIG. 3: Critical boundaries on the mud-rris plane in the NJL 
model (left panel) as compared with the LQCD results with 
Nt = 4 (right panel) . 

In Fig. [3l we display the result for the critical bound- 
ary in our modified NJL model with the finite Matsubara 
frequency summations: N — 15, 20, 50, oo, along with 
the results from the LQCD simulation with Nr =40 
for comparison. We find that the first-order region 
widens as N decreases. This observation qualitatively 
agrees with the arguments in the LQCD simulations 
that the first-order region tends to get bigger than 
its genuine size due to the finite UV-cutoff effects. 
Quantitatively, however, the critical boundaries still look 
far from the LQCD results at Nr = 4. We comment that 
not only the temporal UV-cutoff but also the spatial 
(three-momentum) cutoff leads to similar alteration in 
the critical boundary i.e., we have found that lowering A 
also pushes the critical mass up, though such effects are 
only minor. 

The PNJL model extension It is also worth stud ying 
the critical boundary using the PNJL model BJllJia 
This is because, in view of the results in |22|.|23||. 
the location of the critical point is significantly moved 
toward higher T by the inclusion of the Polyakov loop 
which suppresses unphysical quark excitations. Further- 
more, the PNJL model is quite successful to reproduce 
the LQCD observables for QCD thermodynamics such 
as the quark number susceptibility, the interaction mea- 
sure, the sound velocity, etc. Because the PNJL model 
nicely grasps the essential nature of the QCD phase tran- 
sitions, it would be rather mysterious if only the critical 
boundary on the mud-nT-s shows a huge discrepancy. 

The model is defined by the following Lagrangian 



Cp^jL = Co + Ci + Ce+ll{^,'i>*,T), (9) 

£o = 9 (i^ - 174^4 - n^) g , (10) 
U{^,^*,T) = -feT{54e~"/^$$* 

ln[l - 6$$* + 4($3 + 3) - 3($$*)2] } , (11) 

where U is the effective potential in terms of the Polyakov 
loop. Here <& and $* are the traced Polyakov loop and 
the anti-Polyakov loop which are order parameters for 
deconfinement. They are defined by $ = {l/Nc)trL, 
= (l/iVc)trLt with L = 7'exp[i J^dTA4]. We set 
the parameters a and 6 as a = 664 MeV, b-A^^ — 0.03 
following 0. 
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FIG. 4: Critical boundaries on the m^d-nis plane in the PNJL 
model (left panel) as compared with the LQCD results with 
iV^ = 4 (right panel) . 

Because the Polyakov-loop part is relevant only at 
finite T (as high as comparable with Tc), we can simply 
employ the same parameter sets, as fixed previously 
in the NJL model case with finite TV (see Tab. HI. 
Then, we show the results for the critical boundary 
determined by the same procedure in Fig. |4l Here 
again, we see expanding behavior of the first-order 
region with decreasing TV. Although the PNJL results 
are quantitatively different from the NJL ones, the 
qualitative tendency is the same and, besides, the 
difference is not really substantial. The fact that the 
PNJL model results are such close to the NJL ones is 
in fact non-trivial. The Polyakov loop generally governs 
the thermal excitation of quarks near Tc and particularly 
unphysical quarks below Tc are prohibited by small $ 
and $*. This leads to the consequence that the critical 
point temperature is lifted up to twice as high. From the 
comparison between Figs. [3] and [H gives us confidence 
that the smallncss of the first-order region is not simply 
an NJL- model artifact with unphysical quark excitations. 

Concluding remarks Our model analysis yields in- 
triguing results that the critical masses become larger 
with truncation in the Matsubara frequency summation. 
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This effect is not large enough to give a full account for 
the quantitative difference in the critical boundaries be- 
tween the LQCD and NJL model studies. Nevertheless, 
it raises an interesting possibility that the critical masses 
in the LQCD simulation become smaller if the LQCD 
simulation accommodates larger Nr- This is an impor- 
tant message since the LQCD simulation at Nt = 4, 6, 
and even 8 may have a potential danger that the LQCD 
simulation overestimates the strength of the first-order 
phase transition, and thus the critical point emerges at 
small density due to lattice artifacts. 

An interesting extension of our study is to determine 
the chiral phase transition boundary on the plane with T 
and the baryon (quark) chemical potential ^. It is a com- 
mon folklore that the QCD critical point should appear 
at some Te and ^e- However, it is highly non-trivial to 
determine the order of the phase transition at non-zero /z. 
Recent LQCD simulations at small (i.e. expansion as a 
power of /i/T) disfavor the existence of the critical point 



in the QCD phase diagram 0] . There have been the- 
oretical efforts to understand this LQCD finding in the 
context of the chiral effective model approaches [23 . [25l | . 
However, the errors of the LQCD simulations grow bigger 
at higher /i and, to be worse, the (P)NJL-model analyses 
can be trusted only at a qualitative level so far once /i 
gets large. It is still inconclusive, therefore, whether the 
critical point actually exists or not. Our model analysis 
shows that increasing Nr leads to smaller critical quark 
masses. This result might indicate to further minify the 
possibility for finding the QCD critical point using the 
LQCD simulation unless the discretization errors are un- 
der good theoretical control. 
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